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Abstract 



A system of self-consistent integral equations for the superconducting gap 
is formulated and solved taking account of the instability of the normal phase 
of bilayered cuprates against charge-density waves. The critical parameters 
are calculated as a function of the wave vector, temperature, and doping in- 
dex. It is found that the region in which superconductivity coexists with 
d-type charge-density waves depends strongly on the doping index. The ef- 
fective energy-gap parameter, determined as the interval between the peaks 
of the density of states, can have a local minimum at temperatures T<Tc. 
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It can now be regarded as an established experimental fact that d-type pairing is re- 
alized in most layered cuprates, and in addition the determining component of the order 
parameter of the pseudogap in the normal phase of underdoped samples possesses the same 
d-type symmetry.^'^ At least this conclusion does not raise any objections for YBa2Cu307_5, 
Bi2Sr2CaCu208-5, and other bilayer cuprates. The question of the nature of the pairing and 
the origin of the pseudogap in the normal phase is still open. Specifically, it is still not clear 
whether or not they are interrelated or completely independent. In the present letter we 
show that the empirical dependences of Tc and the pseudogap closure temperature T* on 
the doping index and the dependences of their order parameters on the wave vector can be 
qualitatively understood in a model of a quasi-two-dimensional metal with strong electronic 
correlations, in which case the gap in the elementary excitations spectrum is due to anoma- 
lous averages of the type (^'k'^^'k+q); characteristic for a Peierls instability of the normal 
phase. The problem of the existence of superconductivity under conditions of Peierls insta- 
bility is a classic problem. For ordinary metals it has been discussed in a number of works 
(see, for example, Refs. 3 and 4), where the anomalous averages were due to the Frohlich 
interaction and did not depend on the wave vector. In the present letter the main interaction 
responsible for the formation of the superconducting state is the short-range superexchange 
interaction, and the normal phase exhibits strongly non- Fermi-liquid behavior. The spectral 
weight of the conduction band under study depends nonlinearly on the doping index, and 
it is already half-filled when the number of holes per unit cell of the CU2O4 bilayer is 2/7 
(Ref. 5) instead of 1, as would have happened in the case of normal metals. 

To describe the subsystem of current carriers in bilayer high-Tc superconductors we 
employ a Hamiltonian of the form 

H = Y: e,^r + E C*^'^*^^ + E Jim [2 {SlSrn) - ^] + E 9lmSl5rn, (l) 
1(7 Im l>m l>m 

where ^^(^^^), are Hubbard-type operators corresponding to creation (annihilation) of 
quasiparticles in the bonding singlet-correlated oxygen band,^ rii = -|- "^l^, 61 = ni + 
2^p, and, jim and gim are the superexchange and Coulomb interactions parameters. The 



quasiparticle dispersion relation has the form 

= -Pp(i[2ti(cos k^a + cos kyO) + 4t2 cos k^-a cos kj^a + 2t3(cos 2kj.a + cos 2kya)] — n, (2) 

where is the chemical potential, which under optimal doping, according to photoemission 
data,^ lies 10 meV below the saddle peak of the density of states, Ppd = (2 + 5)/4 is the 
average value of the anticommutators {^^, '^^^^ tunneling taken into account, 5 is 
the number of current carriers per unit cell in the CU2O4 bilaycr, a is the lattice constant, 
and ti = tf'\l + (SiSj) / Pp^. We take the following parameter set, which is consistent with 
the Fermi surface and the temperature behavior of the spin susceptibility of the normal 
phase of YBa2Cu408: tf^ = 70 meV, = 0, and = 5 meV.« Following Ref. 9, we 
give the spin correlation functions of the nearest neighbors phenomenologically as {SiSj) = 
(—0.2 exp(0.6 + S/2)~^ + 0.6)/1.8356, so that the width of the band would approach zero as 
S ^ 0, which is necessary for a correct description of the insulator-metal transition in these 
compounds. 

In the mean-field approximation the elementary excitations spectrum of the Hamiltonian 
(1) is determined by the equation 

-ev-E 
-Gt 



det 



e-k — E Gk+Q 

Ck £k+Q — E 

At 



TI* 

•^k+Q 



C/k 

Ak+Q 

-Gl 



Ut 



^k+Q 



'k+Q 
-£k+Q — E 



0. 



(3) 



The order parameter of the superconducting transition (ST) is given by the expression 
1 



E[2j(ki - k) - g{k, - k) + P^,aB{ku k, - k)]{^\^'^ 



(4) 



where j(q), 5'(q), and i?(ki,q) are the Fourier transforms of the superexchange. Coulomb, 
and Frohlich interactions potentials, respectively. Specifically, 

i(q) = io(l - 65^) ■ {cosq^a + cosg^^a), 

where jo = 125 meV for YBaCuO compounds^" and 



5'(q) = exp(-7(5) (cos q^a cos qyO) , 



where go — 315 meV. Among all possible phonon modes, the buckling modes make the most 
important contribution to pairing. In this case 

9 1 + i fcos (/a-a + COS o^al , , , , 

B{k, q) = 2Eg^cUq ^ „ 7 ^-\^e{huD - kk - e^^UMhuD - kkl), (5) 

[nujci) — (,£k — £k+qj 

tiwn is the Debye frequency of the order of 45 meV, Q{x) is the theta function, and Bq ~ 
35 meV. In the range of doping indices S of interest to us, Eq. (4) has a solution of the form 
Ak = Ao(cosA;a;a — cos kyu). 

As was pointed out in Ref. 5, the CDW parameter Gk is a sum of two components with 
a different dependence on the wave vector. The component associated with the short-range 
potentials 

Gt: = Eb-(ki - k) + g{k, - k)](*^%IjQ), (6) 
possesses d-type symmetry — iGo{coskxa — cos kyo), while the phonon part 

possesses s-type symmetry. It is relatively small and is neglected below. The order pa- 
rameter Uk entering in Eq. (3) is related with thermodynamic average (^k^^I^k-q)' which 
appears naturally in the equations of motion in the presence of superconductivity and Peierls 
instability. In the general case it is determined by the expression 

= -W-^ Eb'(ki - k) + j(ki - k - Q) - g{k, - k) + P^M^i, ki - 

^pd^^ k 

(8) 

We investigate below the approximately optimal doping regime, where Q ~ (tt, tt). One can 
see from Eq. (8) that is comparatively small and can possess only s symmetry, since for 
Q — (tt, tt) the sum _7(ki — k) -|- _7(ki — k — Q) vanishes. We note that for C/k = Eq. (3) 
leads to the spectrum 



£^1,2 = \ \{el + 4h-q) + AkAJ^ + G^Gl ± \eI^ (9) 



where 



(10) 



and £'4 = — and — —E2. 

In the general case we have the following system of integral equations: 
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where 



TVa = ^^{el - 4+q) + AkG^Gk + G*) + t/k(Akt/i: - A* t/k) + 2£k+QG;;t/k, 
A^G = ^Ckl^k + £k+Q)' + AkA* (Gk + Gl) + £k+QA* t/k - £kAkC/^, 
A^C7 = ^C/k(£k - £k+Q)' - Ak(Ak[/i; - A* C/k) - £kAkG* + Sk+qAkGk. 



(12) 



The system of equations (11)-(13) was solved numerically. As expected, the parameter 
Uk is neghgibly small. The computed phase diagram of the S dependence of the critical tem- 
peratures for Ao(Tc) and Go{T*) at Q = {^tt, ^n) is close to the corresponding diagrams 
proposed in experimental works (see, for example, Ref. 11). In Fig. 1 three different doping 
cases are presented as an example of the temperature dependences found for the amplitudes 
of the order parameters Ao(T) and Go{T): ai, 0:2 and 0:3, where the parameter a is given 
by the ratio a{S) — {S — Smin)/{^opt ~ ^min)- Here Smin is the minimum value of the doping 
corresponding to the appearance of a superconducting state. In all cases, for 5 < 5opt — 
0.315 (i.e., the existence of a CDW) the temperature dependence of the superconducting 
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order parameter differs considerably from the BCS case. In the case of a strongly overdoped 
regime the ratio 4Ao/fcsT(7 — 4.2; for 5opt it equals 4.5; it then increases, as seen in Fig. 1, 
and reaches values of the order of 10 for 5 ^ 5min (i-e., in the strongly underdoped region). 
This kind of variation of AAo/ksTc has already been noted in a number of experimental 
works. The computed temperature dependences of the pseudogap Gq exhibit nonmono- 
tonic behavior as temperature decreases; this has not been noticed previously. It would 
be interesting to check this prediction of the theory experimentally. Thus, if the effective 
gap parameter Ae//(r), determined as the energy interval between peaks in the density of 
states, is measured in scanning tunneling spectroscopy experiments (see, for example, Ref. 
12), then on account of the nonmonotonic behavior of Gq{T) anomahes should be observed 
in the plot for Aeff{T): 1) In type ai samples Ae//(T) will decrease (!) as T ^ 0; 2) for 
q; ~ as a local minimum should be observed in the plot of Ae//(T). Examples of such plots 
calculated on the basis of expressions (9) and (10) are presented in Fig. 2. In our opinion 
the indicated anomahes can be seen from the experimental data.^^'^^ These results, however, 
require a more detailed experimental check on samples with different doping indices. 

This work was supported in part by the Russian State Science and Technology Program 
"Superconductivity", Project No. 98014. 
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FIGURES 
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FIG. 1. Computed temperature dependences of the amplitudes of the order parameters Aq 
and Gq. The curves corresponding to the superconducting gap parameter Aq start at the point 
T/Tc = 1. The relative-doping parameter a is defined in the text. 
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FIG. 2. Predicted examples of nonmonotonic temperature variations of the effective gap Ag// 

between peaks of the density of states of the energy spectrum. 
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